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We study in a Brill-Hartle type of approximation the back reaction of a superposition of linear 
gravitational waves on its own mean gravitational field up to second order in the wave amplitudes. 
The background field is taken as a spatially flat Einstein-de Sitter geometry. In order to follow 
inflationary scenarios, the wavelengths are allowed to exceed the temporary Hubble distance. As 
in optical coherence theory, the wave amplitudes are considered as random variables, which form 
a homogeneous and isotropic stochastic process, sharing the symmetries of the background metric. 
A segregation of the field equations into equations for the wave amplitudes and equations for the 
background field is performed by averaging the field equations and interpreting the averaging process 
as a stochastic (ensemble) average. The spectral densities satisfy a system of ordinary differential 
equations. The effective stress-energy tensor for the random gravity waves is calculated in terms 
of correlation functions and covers subhorizon as well as superhorizon modes, where superhorizon 
modes give in many cases negative contributions to energy density and pressure. We discuss solutions 
of the second-order equations including pure gravitational radiation universes. 
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I. INTRODUCTION 

Many classical aspects of gravitational waves are still 
poorly understood. Most of them are connected with 
the fact that the full theory of gravitational radiation 
as following from general relativity is extremely nonlin- 
ear. Surprises can be expected, if the nonlinear regime 
becomes more deeply explored. Pure numerical meth- 
ods acting on the full field equations [jl) are important 
but cannot be the only key to the nonlinear regime. 
They should be supplemented by suitable approximation 
methods, which allow an analytical or semi-analytical 
approach. Already in 1964 Brill and Hartle B pro- 
posed a scheme, which takes the back reaction of lin- 
ear gravitational waves on the background metric into 
account. Usually, the Brill-Hartle method (see also [pi, 
B) is considered as a high-frequency approximation for 
gravitational radiation. For cosmological application the 
full spectrum of gravitational radiation including low- 
frequency (superhorizon) modes must be studied, if its 
generation and propagation through inflationary stages 
is considered. We discuss in this article a different inter- 
pretation of the Brill-Hartle approach, which allows to 
treat also low frequencies. 

The Brill-Hartle method requires an average over small 
scale ripples of the geometry. How averages can be for- 
mulated consistently in general relativity, is a basic and 
still not completely solved problem [pi. Our assumption 
is to interprete the perturbations of the geometry as ran- 
dom functions as in optical coherence theory B or in the 
Monin-Yaglom approach to statistical fluid mechanics [J7) . 
The averages of the Brill-Hartle method are then taken 
as ensemble averages, not as averages over space-time 
regions. How stochastic averaging and a random geome- 



try can be reconciled with general relativity in a rigorous 
way is also an open question. The problem is not touched 
upon here, we adopt instead a field theoretical approach 
in the sense of the Monin-Yaglom treatment. 

Along this line we discuss the response of the expan- 
sion rate and of the wave amplitudes to the mean grav- 
itational field produced by the waves. We assume that 
the stochastic process for the wave amplitudes is homoge- 
neous and isotropic. It is not necessary to suppose it to 
be Gaussian, since only two-point correlation functions 
are involved at the approximation level considered here. 
The background is kept as simple as possible, a flat 3- 
geometry with scale factor a(rj) and conformal time w is 
assumed. The treatment is completely classical, we do 
not consider back reaction effects from quantum grav- 
ity ||, semiclassical gravity |9j or an Einstein-Langevin 
equation |1(|. With methods similar to those discussed 
here, the subject was treated in |0| and |P^[ . 

The paper is organized as follows. In section II the 
segregation of the field equations into equations for the 
wave amplitudes and equations for the background field 
is discussed and, following earlier treatments Q, Q, 
the effective stress-energy tensor for the random gravity 
waves is calculated in terms of correlation functions for 
the wave amplitudes. Depending on the wave spectrum, 
major contributions to this tensor may come from waves 
with wavelengths exceeding the horizon distance. Sec- 
tion III is concerned with gauge problems. The rest of 
the paper treats solutions of the averaged field equations 
with various assumptions for the spectrum of gravita- 
tional waves and for the presence of matter fields. The 
aim is to discuss effects of back reaction in general, we do 
not consider the origin of the waves and will also not take 
observational constraints into account in this article. In 



section IV we consider the evolution of the gravitational 
wave amplitudes for different frequency regimes in cosmo- 
logical models, using a number of simplifications. Pure 
gravitational wave models ("geons") are the subject of 
Section V. Its discussion is complicated by the necessity 
to include a non-linear back reaction term in the wave 
equation. In the final section VI possible improvements 
of the approach are shortly discussed. 



II. AVERAGING THE FIELD EQUATIONS. 

We take an Einstein-de Sitter model as background 
metric and add tensor perturbations, writing^ (i,k = 
1,2,3) 



ds 2 



-a(rj) 2 dri 2 + g ik dx l dx k 



with 



9ik — fl &ik + hik, hu — 0, h ik _ k — 0. 



(1) 



(2) 



It is appropriate to write the field equations to this metric 
in the (3+1) form, assuming a perfect fluid as matter in 
the background metric: 

r(3) +k 2 - K lk K lk = 16TrGp m , (3) 



9 Kl K, 



hill 



K, 



0, 



-K'fk/a + i4 3) + KK lk - 2K a K k l = 
4nGgik(p m — Pm)- 



(4) 



(5) 



The covariant derivatives are denoted by a prime and 
taken with respect to gi k , and ' = -4-. Indices are moved 

with the three-metric gi k . R ik and R^ are the three- 
dimensional Ricci tensor and scalar for g^ The extrinsic 
curvature Ki k is given by 



K ik = ~Ya 9 ' ik = ~ a ' Sik ~2a h ' ik ' 
Up to second order we have 

g ik = S ik /a 2 - h ik /a 4 + h u h k i/a 6 . 



(6) 



(7) 



With the decompositions (0) and (0) one can write 
down the field equations explicitly, including all terms 
up to the second order. Simplification results from the 
gauge restrictions for hik, which are taken into account. 

(3) 

The Ricci tensor R\ k is given by 



»( 3 ) 1 h 

Rik ~ ~^2 hlk > 11 

T~ T iLlm\" j ik.lm > i^lmAk ^iLkm ^kLim) 

+ -T—lhl m ,ihl m! k + -^—[(hil im hkl,m — hu jm hkm,l)i (8) 

and the field equations read explicitly 
13 1 

+ —^(hklhkl,mm + -rhkLrnhkLm ~ -zh k i m h km i) 

a° 4 2 

a' a' 2 1 
f h k ih' kl + 3—h kl h kl - —^h' kl h' kl = l&irGp m , (9) 



4a 6 



- ^-Qhkihki,i gh k ihik,i + ■^-ghkih ik<l 

-^h kl ti kll + —h' kl h k i,i - -^hki,ih' kl = 0, (10) 



{— + ~^> dik + ~^2 h ih - ~^ h ik + 2 -X«* 



a' 2 1 , , a! , a! , 

-2-jhuhki - 2-±h u h kl + -rh u hu + —gh ki hu 



a" 



a' 
2a~ 5 ' 



a" 



(3) 



— -^ZK^ ik ^lm^lm ^ QSikhi m hi m + i? ife 



4nGgik(ftm-Pm)- (11) 



The Brill-Hartle method usually starts with the assump- 
tion that the space-time variation of the small ripples 
h^ (of order h/X, where A is a typical radiation wave- 
length) is much larger than the variation of the back- 
ground metric (of order 1/L). Thus terms in the Ricci 
tensor which are bilinear in h^p (of order vr) are com- 
parable to terms involving the background metric (of or- 
der -jtj), if h rs y. This latter condition excludes low- 
frequency waves A > L, since h must be sufficiently small 
to allow a second-order approximation. Terms linear in 
hjiv,p<r are much larger and should therefore vanish sepa- 
rately, giving rise to the linear wave equation for h^ v . 

We do not follow this bookkeeping ( it was criticized 
in Jig]). Instead, the functions hik are interpreted as 
random functions. Performing a stochastic average of 
([|[nj[ll]) removes terms linear in h ik , but keeps terms 
bilinear in hik and in the derivatives of hik ■ This allows a 
segregation of the field equations without restrictions for 
the wavelengths. The result of the stochastic average can 
formally be written as (Eqn.(|l0|) reduces to an identity) 



@ 



We follow the conventions by Thorne, Misner and Wheeler 



a 
a 



3-f = 87rG(p m + Pg ), (12) 

a' 2 
+ — r = 4:nGa 2 (p m -p m + p g -p g ), (13) 



where p g and p g are the averages over nonlinear terms. It 
is convenient to interprete p g ,p g as the effective density 



and pressure of the gravitational radiation field. Sub- 
tracting the averaged field equation ( p~3| ) (multiplied with 
5ik) from ( pd| ) and neglecting higher-order terms of the 
type (A) — A, as is usually done in a self-consistent field 
approximation, one obtains a modified wave equation for 
the amplitudes hik 



ik 



Ah lk - 2—h lk + 2h ik ( 



a' 2 a 



with 



b = 4TT(Ga 2 (p g ~ Pg ). 



+ 6) = (14) 



(15) 



This equation differs from the usual form by a time- 
dependent term b representing a back reaction of the 
energy and pressure of the waves on their propagation 
(a factor of £ was introduced in front of b in order to 
switch off the back reaction term for comparison pur- 
poses) . Since p g and p g depend on solutions of the wave 
equation, some nonlinearity is thus introduced, b re- 
sults from second-order bilinear terms and is therefore 
neglected in linear treatments of cosmic wave propaga- 
tion, but should be kept in the spirit of our approach. In 
situations where the waves do not appreciably influence 
the scale factor evolution, the back reaction term can 
be neglected. Also its influence is usually small in the 
high-frequency approximation, when the wavelengths are 
small compared to the Hubble distance. Note also that 
the back reaction term is exactly zero, if the wave back- 
ground has the Zeldovich equation of state p p — p g . On 
the other hand, if low-frequency radiation contributes ap- 
preciably to the average density and pressure and hence 
to the scale factor evolution, the modified wave equation 



must be considered in general. (A previous paper ( |17|) 
on the same subject was based on the linear wave equa- 
tion with b = 0. We shall continue its use in Sec. IV 
in order to compare later the results with those of the 
general case b ^ 0). 

Some authors use different definitions of the wave am- 
plitudes by applying factors of a on them. We have de- 
fined the spatial components hik as perturbations to the 
three-metric a 2 Sik- Equivalent to (111) are 



(h lk /a 2 ) -A(h tk /a 2 ) + 2-(h lk /a 2 ) + 2bh lk /a 2 = 

a 



(16) 



and 



{hv, /a) -A(h lk /a) + (2b ){h ik /a) = 0, (17) 

a 

which are sometimes easier to use. It was fre- 
quently noted that gravitational wave perturbations in 
a Friedman-Robertson- Walker (FRW) universe may be 
described as a pair of massless minimally coupled scalar 
fields in the same background space-time (see, e.g., Ford 



and Parker in ( M). If back reaction is important, this 
correspondence is lost in general, but for a de Sitter scale 
factor a ~ 1/rj the modified wave equation has the form 
(□ + £R + m 2 )(<$>/a) = with £ = | - ±m 2 aV, which 
characterizes non-minimally coupled and possibly mas- 
sive scalar particles. 

To calculate the components of the gravitational stress- 
energy tensor, we assume that the random process rep- 
resented by hik is homogeneous and isotropic. The cor- 
relation functions which enter (0,0) are related to cer- 
tain spectral densities. To see this, we represent hik as 
stochastic Fourier integral (see B for a detailed treat- 
ment of the spectral representation of random processes) : 

hik(x,rf) = j lk (k,r))e lWx dk + conj.compl. (18) 

From (|l4|), the Fourier amplitudes satisfy the ordinary 
differential equation 

7a " 2 77- fc + 7^(2^ -2^- + k 2 + 26) = 0. (19) 

Amplitudes of the correlation functions may now be writ- 
ten as frequency integrals over spectral densities: 

(^ fc (x, r))hi m (x, r])) = / ((jikim + "f* k i m ))dk, (20) 



whereQ §} 



(7ifc(k,r?)7f m (k,r?)) = <5(k - k) 7ifcim , 



(21) 

(22) 



It is easy to see that the spectral densities jikim satisfy 
the symmetry relations 



* 

liklm — flmik i ^iklm — ^kilm j 
likll = 0, Jiklmk" 1 = 0. 



(23) 
(24) 



The general solution of the algebraic constraints ( 23[ 24) 
for the spectral densities contains four complex functions 
of k and rj, which describe polarized background radi- 
ation in general. The stochastic background of gravi- 
tational waves expected from pre-galactic stages of the 
Universe could be polarized due to strong anisotropics 
expected at Planck time. Here we confine the discussion 
to unpolarized radiation, represented by a single (real) 
spectral function a(k,rj). Then "/ikim can be written as: 



liklm = Ct{k,T])Siklm, 

Siklm = SuSkm + SimSkl ~ 5i k 8l m , 



(25) 
(26) 



' We refer to equal space - equal time correlators of the form 
(poh as "correlation functions". 



(27) 



where 

Oik = uik fcikk j «* 7 & = K K . 

The transversal <5-symbol satisfies 

4 = 2, 4fc fe = 0, 44; =4- (28) 

To calculate the averaged stress-energy tensor, we need 
also the spectral densities 



(likllm) = P s ( k ~ k)5 ik im, 

{likl'Cn) = 7^(k - k)Siklm, 



(29) 
(30) 



where /3 is real and 7 is complex in general (the correla- 
tion functions a, fi and 7 were introduced in J14|). Simi- 
lar to (pO) , the correlation functions which are needed to 
calculate the effective wave stress-energy tensor can be 
written as 



\'lik,r'llm,s) — ^ I ^Oiklm^r^s^'^^ 
\llik'l'lm,rs) — \i^ik,ri^lm^s) 5 
( h ik h 'lm) = 2 / PSiklmdk, 



(h ik ti lm ) = 2 / 5R(7)4 im dk, 
{hikhim.r) — 2 / a/c r 4-; m dk, 



(hikh' lmr ) = 2 ^(rfSikimkrdk 

{hikhim.r) — ~\hikhi mr ). 

Since a, /3, 7 are functions of k (and 77) only, the angular 
integrations can be performed easily. This gives with 

Piklm = 44m + SimSkl- 

(hikhim) = -——{ipikim - 2S ik 6 lm ) / ak 2 dk, 

(hik,rhlm.s) — 77— ■(— 10&tfe4n4s + HPi/rim4s 

105 

^Pilrs^kra ^Plkrs^ira 
dPkmrsVil ^Pimrs^kl 



(hikhi. 



(h'ikh'i m ) = -r--(3pifej m - 24<5; m ) / /3£; 2 dfc, 

(h lk ti lm ) = -— (3p jfc z m - 24<W / K(7)fc 2 rfA;, 
(hikhim^) = 0, 

(hikhlm,r) = 0, 

{hikh' lmr ) = 0. 

Since only the real part of 7 enters the averaged ex- 
pressions, we denote subsequently this real part by 7. 



Some unusual results can be expected when one deals 
with stochastic averages over nonlinear quantities. For 
the averaged three-dimensional Ricci tensor we obtain, 
using the just derived relations, 



<^4 3) > = |^4 J ak A dk. 



(31) 



Thus, whereas the three-dimensional background met- 
ric is flat, a random superposition of gravitational radi- 
ation to this background produces an averaged 3-Ricci 
tensor with positive curvature, if curvatures are defined 
as eigenvalues of the Ricci tensor. On the other hand, 
the stochastic average of the 3-curvature scalar R^ = 
g lk R ik is negative: 



(ij(3)) = _??E f a kUk. 
a b J 



(32) 



The bilinear terms which enter R^ and produce a 
nonzero stochastic average are different from those en- 
tering Rrj^, thus the result (Ri k )(g %k ) ^ {R lk g lk ) docs 
not come as a surprise. But it makes it hard to inter- 
prete stochastic averages geometrically. 

Energy density and pressure of the gravitational waves 
are given by similar integrals 



Pg = 



2Ga 6 



Pg 



l-g [ dkk 2 {k 2 a + /3 + 4 7 - - 12cA-), (33) 

\— [ dkk 2 (7k 2 a-5[3 + 20j—-20a'^-). (34) 
xa b J a cr 



6Ga 6 



The deviation from the HF equation of state follows as 



p g - ip g 



1 



Ga 6 



,'2 



dkk 2 {-ik 2 a + 3/3 - 87— + 4a-5-), 



and the back reaction function b is given by 



8ttC 
3a 1 



dkk 2 (-k 2 a + 2/3-2~f- 



4a- 



(35) 



(36) 



The time evolution of the spectral functions a, /3 and 7 
is obtained from ©. Differentiating @,(§9]),(|^) with 
respect to 77 and using (|19|) leads to a system of coupled 
differential equations: 



27, 



.'- 



= 4-/3 + 2 7 (2- - 2\ - 2b - k 2 ), 



a 



a a 

11 



7 



a , a 

2— 7 + /3 + a(2 — 

a a 



26 -r 



(37) 
(38) 

(39) 



With given initial values the time evolution of all spec- 
tral densities and hence correlation functions follows from 
this system, if the scale factor is known. The system for 
a, /3, 7 is actually nonlinear due to the presence of the b 



term. It can be transformed into a single nonlinear dif- 
ferential equation for a alone. Substituting a'/2 for 7 in 
the last two equations and introducing a new function e 
instead of by means of 



11 



e = /3- 



a 

4a' 



one obtains 



, , a a. 

e'+ 4— )e = 0. 

a a 



This is integrated to 



e (k)c 



(40) 



(41) 



(42) 



Then the basic spectral function a satisfies the nonlinear 
differential equation 



over the spectral density f(k,rf), denoted as "moments'' 
subsequently: 



/o = / dk k 



eo(fc) 



f(k,r,Y h 



dk k 



2 nM) 



h = / dk k 2 f(k, V ), h= dk fc 4 /(M). (48) 



In the density and pressure equation, / and /1 appear 
only in the combination 



.94 = /o + A/4. 
We also note the useful relations 

ffi = -/i+(--2&)& 

a 

/£' = 2.g 4 + 2(— - 26)/ 2 - 2/4 



(49) 

(50) 
(51) 



2aa" - a' 2 - 4e a 4 - 4-aa' - 4a 2 (2 2— - 26 - fc 2 



= 0. 



(43) 



If we define wave amplitudes with a different power of a, 
for instance as in (|l7|), this equation can be simplified. 
With a = fa 2 , 

VI" - f' 2 + 4/ 2 (k 2 + 26 - — ) - 4e = 0. (44) 

a 

This equation is also nonlinear in /, but its nonlinear- 
ity arises because we deal with the spectral density of an 
expression which is quadratic in the random wave ampli- 
tudes hik ■ The solutions of (M) are related to those of a 
differential equation, which is linear in the wave ampli- 
tudes. Let hi, hi be a pair of real solutions to 



such that 



h" + h{k 2 + 26 - a" /a) = 



h\h' 2 — h%hi 



eo, 



(45) 



(46) 



where eo is a time-independent function of k only. Then 
/ = h\ + h\ is a solution of (|44| ) . Not surprisingly, the lin- 
ear differential equation for h is the Fourier transform of 
the wave equation ( |l7|) for a single realization of the ran- 
dom process. We also note that the solutions of (|44| ) can 
be reduced to solutions fhom{k,rj) of the homogeneous 
part of ( f44|) (eo = 0): Any solution of the full equation 
( [44] ) can be written in terms of a suitable function fhom 
as 



/ = fhom +£o(k)f ho 



di] 



fi 



home 



(47) 



It is convenient to write the expressions for p and p and 
for other stochastic averages in terms of four frequency- 
independent, but in general time-dependent, integrals 



which follow from differentiating f\ once and fi twice 
and using the differential equation for /. 

For a general scale factor, energy density and pressure 
may be rewritten as 



P 9 = ^(U + 54 + 3-/^-7^/2), (52) 

zGa 4 a a z 



3p ^2ij (7/4 - 5 - 94+5 7 /2 ~ 5 ^ /2) 

The back reaction function 6 is given by 

*=It(-/4 + 2.94 + -f 2 -^h\ 
6a A a a z 



(53) 



(54) 



thus the relations (pO) , (pll) for the moments become non- 
linear in general. 

For completeness, we give the ensemble averages of 
some other geometrical quantities. The mean value of 
the four-dimensional Ricci scalar can be written as 

(R) = K + ^( 3 /4 - 3 34 + ^h + -&\ (55) 
a° a* a a 

where the first term is the background contribution. The 
second term may also be written as 8irG(3p g — p g ). Since 
R = in vacuum, also the stochastic average of R must 
vanish. Then ( J55[ ) is equivalent to a combination of (|1J) 
and (O) (with p m = p m = 0), as one can check easily. 
There are only two independent and in general nonzero 
components of the averaged Riemann tensor: 



(i?oioi) = a' 2 - aa" + ^(g 4 + ^/ 2 - -/^), (56) 
3 a z a 



\Tx 



a' 2 . a' 



(Ri2i 2 ) = a' 2 + —(3/ 4 - ff 4 - — h - -fit)- (57) 
3 a z a 

The one independent component of the averaged Weyl 
tensor 



(O)10l) - (Cl212) = ^(/4 - ^h + f&)- (58) 
d a z la 

contains, of course, no background term. All curvature 
components result from averaging nonlinear terms, they 
have no direct connection to the first order perturbations 
of the curvature tensor, whose average is zero. It is im- 
mediately seen that {g a ^){C a ij,0v) differs from zero for 
some components, while, of course, the original expres- 
sion g a ^C ai _ t p L , and hence its average vanishes. Averaging 
the four local curvature invariants connected with the 
Weyl tensor is more complicated: Apart from the cor- 
relation functions listed above further expressions such 



h ik (x i ) = h ik (x i )-a 2 (e k + &) 



(61) 



as (h 



ij,kl "mn,rsj 



must be calculated, which introduces 



higher-order moments such as fa = Jdk k 6 f(k,rj). 

We summarize the main results of this section. In or- 
der to describe the back reaction of gravitational radi- 
ation propagating in a flat Friedman universe upon the 
scale factor, we have replaced the Einstein field equations 
by a system of equations for averaged quantities. This 
system consists (i) of the Friedman equations (|l^,(|l3|) 
for a cosmic fluid with density p m and pressure p m , to- 
gether with the corresponding quantities p g and p g for 
gravitational radiation given by (|5^) , ( p3| ) , and (ii) of the 
differential equation (H) for the spectral density f(k,r]). 
The latter equation may be replaced partially by differ- 
ential relations for the moments of /. The aim is to solve 
these equations for the scale factor and for / simultane- 
ously. One has to keep in mind that the results are only 
reliable if terms higher than second order in the original 
field equations can safely be neglected. 



III. GAUGE TRANSFORMATIONS 

We have chosen to work in a fixed gauge. As re- 
cently noted by Unruh llq], there is no basic difference 
between a gauge invariant approach, where the equa- 
tions are written in terms of "gauge-independent" vari- 
ables (see, e.g., [Oj), and an approach where the gauges 
are fixed in some way. It is only important to fix the 
gauges completely, leaving no room for residual gauge 
freedom, which could be mistaken as physical degree of 
freedom. We therefore ask for the coordinate transfor- 
mations x^ = x^ + £ M preserving the eight constraints 
^oo = 0, hoi = 0, ha = 0, hik.k = 0. To linear order in 
£ M the perturbations transform as 

V(»") = M*") - &? $}& ~ fl$&- (59) 

A simple calculation shows that for a general scale factor 
the transformations we are looking for must satisfy 

e° = o, e = e{x k ), e, 



0, £ 



u- 



0. 



(60) 



The functions £ M therefore do not necessarily form a 
Killing field, and the first-order perturbations are gauge- 
dependent in general. The spatial components transform 

as 



It is easy to see that the first-order wave equation ( fI4| ) is 
invariant with respect to (pTI).The transformation of the 
correlation functions is derived assuming that the vector 
£* is also a random process with zero mean and indepen- 
dent of the wave process. We outline the procedure. The 
random process £* must be homogeneous and isotropic 
as hik, so its correlation function can be written 



(ee 



dk(Fi(fc) 



k l k h 



F 2 {k)5 tl 



(62) 



£\ = is translated into Fi(k) = —F 2 (k) in the spectral 
region. It is more difficult to write the condition ^ l kk = 
directly as condition for the spectral density F\, so 
we first proceed without this constraint. With ( pl|) one 
forms (hikhim), this gives 

J 2 = f 2 + — I dkk i F 2 (k). (63) 

Similarly, the reduction of (h ik h lm ) yields 

g 4 = 94 + ^-Jdkk 4 F 2 (k). (64) 

Writing out (hik, r him jS ) shows that the average of the 
^-depending terms does not have the symmetries of the 
other terms. This reflects the fact that the condition 
£ l fcfc = was not taken into account. We conclude 
that the spectral density satisfies the integral relation 
J dkk 6 F 2 (k) = 0, and that f 4 transforms as 



h = h 



(65) 



It is easily seen that the density p g and pressure p g as 
defined by (]35 ),(p4] ) are invariant against the transfor- 
mations (|63])-(35|), and this holds also for the averaged 
components of the Weyl tensor and for the Ricci scalar. 
As is well known, the effective stress-energy tensor is 
not gauge invariant with respect to general coordinate 
transformations, i.e., transformations which violate the 
constraints. However, as shown by Abramo, Branden- 
berger and Mukhanov ( JO] , see also Q , pj), the gauges 
change at the same time the background geometry to sec- 
ond order, and these changes just compensate the change 
of the energy-stress tensor. 



IV. SOLUTIONS INCLUDING MATTER FIELDS 

We now try to find solutions of the equations derived 
in section II, assuming that apart from gravitational ra- 
diation other forms of matter are present, which domi- 
nate dynamically during most stages of the cosmic evo- 
lution, if not always. Dominance of matter means that 
the 6— term in the wave equation can be neglected. 



A. The high-frequency regime 

The high-frequency (HF) regime is defined by the as- 
sumption that all wavelengths are small compared with 
the temporary Hubble distance. Then the time deriva- 
tives which occur in (|37|-p9|) are small: Let A be one 
of the quantities a, (3, 7, then A' is of the order A/T, 
where T is a Hubble time, T ss a/a' , and this is much 
smaller then a multiplication with the wave frequency k: 
A/T -C Ak. The terms (3 and k 2 a in (B3) are therefore 
much larger than other terms and must cancel in a HF 
approximation: 



/3 = k 2 a. 



(66) 



Using this relation and neglecting time derivatives in (|3£ 
p4|), one obtains 



Pg = 3Pg = 



1 



Ga 6 



dkk a. 



(67) 



equivalent to the equation of state for a gas of noninter- 
acting massless particles. A similar cancellation of terms 
must occur in (|43|). This requires £q > 0, gives 



/eo<2 /k 



(68) 



and identifies the function eo(fc) as closely related to the 
time-independent spectral function in the HF regime. 
The energy density in this regime, 



1 



Ga 4 



dkk -Jcq, 



(69) 



shows the typical 1/a 4 dependence on the scale factor 
for background radiation. Apart from an overall redshift 
factor, the gravitational wave spectrum does not change 
in time. 

The expression ( p7j ) for the effective gravitational en- 
ergy density and pressure of high-frequency gravitational 
radiation holds also in the case that this radiation con- 
tributes appreciably to the background geometry, beside 
other forms of matter. The result is only a changed time 
dependence of the scale factor a. 



B. Stress tensor evolution in a fixed background 

Considering now the full spectrum including both low- 
frequency (LF) and HF modes, one must return to the 
solution of the complete equation (W4) . The change of the 
spectrum during expansion depends on the wavelength. 
The solutions / of (B3) for power law scale factors can 
be represented by Hankel functions of second kind, if we 
neglect the b— term assuming that gravitational waves 
do not contribute appreciably to the scale factor. We 
prefer explicit expressions, since this allows to use Fourier 
methods. For the matter universe (a ~ rj 2 ), the radiation 



universe (a ~ 77) and the de Sitter universe (a ~ 1/rf) as 
background geometry the general solution can then be 
written 

a /a 2 = f = 2npp* + (I + im)p 2 + (I - im)p* 2 . (70) 

/, m, n are three functions of k and connected with eo by 

e = 4k 2 (n 2 -l 2 -m 2 ). (71) 

p(x) is a complex function of x = krj, given by 

p(x) = (1 + is/x)e ix (72) 

with s = 1 in the matter and de Sitter universe and s = 
in the radiation universe. (In spite of the different scale 
factors, the spectral density / in the de Sitter cosmos 
has the same form as in the matter universe, since a" /a 
is the same. This coincidence holds only for the spectral 
density / and is destroyed, if one returns to a). For large 
x, f reaches the same asymptotic form in all models. 
This asymptotic form coincides with the HF expression 
/ = y/eo/fc considered previously, if the ratios — and — 
tend to zero for large k. 

How do the different wavelengths contribute to the ef- 
fective stress-energy tensor? It is usually assumed that 
modes with wavelengths larger than the temporary Hub- 
ble radius can be discarded, since they are unobservable 
and not true waves Efj| , have no appreciable influence on 
the local wave energy density pl|, |22f|, [J23[ , or appear 
locally as gauge transformation |24j| , |18||. The relations 
(p2p and (p3) allow a definite answer. For a radiation 
universe a ~ 77 one obtains: 

7 7 3 

a A Gp g = 2n A -n 2 - jj-^Vkfa) + ^Z^M, ^ 

5 5 5 

3a 4 G Pg = 2n 4 - -^n 2 - l^Mv) + T^^iv) 

-\^{jl) (74) 

with the time function 

^ 2 {v) = 2 / k 2 (l(k)cos{2kri) ~ m(k)sin(2kr]))dk, (75) 
Jo 

we also use the notation ni — L n(k)k l dk for the time- 
independent moments of n(k). The spectral functions 
l,m,n are constrained by the condition that the fre- 
quency integrals (including the moments n 2 and 714) con- 
verge. To justify the assumption of a pre-determined 
scale factor, the amplitudes must be small enough that 
their contribution to the background geometry is negligi- 
ble. If the function ip 2 (77) with its first time derivative is 
bounded for large r\ and ip 2 tends to zero for large 77, the 
relation p g = p g /3 between pressure and density follows 
asymptotically for large 77, as seen from (|73|) , (f74[) . For 
small times however, p g and p g show strong deviations 



from this relation. In particular, p g as well as p g become 
negative for sufficiently small n. In the case I — m = 
(and therefore ip2 = 0) this is immediately seen from ( 173 ) 
and (um . p g and p g are also negative for small n in more 
general spectra. 

The main reason is that the spectrum for wavelengths 
exceeding the local horizon scale at n (corresponding to 
frequencies k < 1/n) gives a negative contribution to the 
total values of energy density and pressure. If one goes 
back in time, the superhorizon modes ultimately domi- 
nate the spectrum (we always assume that the spectral 
density in the ultraviolet region drops sufficiently fast to 
ensure convergence for k — > oo. Ultraviolet divergence is 
a problem of quantum theory and does not concern us 
here) . 

Since the HF part of the spectrum satisfies a p g a 4 = 
const-law and is therefore not very interesting, we con- 
centrate on the LF tail, which is essential for back reac- 
tion effects. 



C. Back-reaction solutions in the regular 
low-frequency limit 

In the long-wavelength or LF limit JnOwe assume 
that the spectral density / can be expanded in powers of 
k around k = 0. This assumption excludes a singularity 
at k = 0, which in some cases can be interpreted as lead- 
ing to a finite infrared contribution to p and p (Section 
V). Writing down (44) explicitly requires to calculate the 
back reaction term b and hence the integrals in (Hq). In 
the sense of our approximation, this is trivial if the in- 
tegration is assumed to extend to a maximal (but still 
small) k = k\. Expanding also eo, (J44]) leads to 



2qq"-q n 



327rfc|C , 1,2 2 , 1 , 2jL ,a' 
6 5 6 a 



Aqq'- 
a 



(76) 



as lowest-order equation of a hierarchy of equations 
for the expansion coefficients (we have written q for 
limfc^o f/a 2 and again eo for lim/j^ofo)- Only this 
lowest-order approximation will be discussed. Since the 
solution of ([44 ) can be written in the form (|47|) , it is ap- 
propriate to treat the case eo = 0. Here q — const is an 
obvious solution of (|76|), provided we neglect the term 
proportional to k\ in the bracket, which is small com- 
pared with the other terms in a LF approximation. The 
effective energy density and pressure of the extreme LF 
background follows as 



Pa 



3 _ qk\ 



(77) 



Note the a 2 decay compared to the a 4 decay of the 
energy density of the high-frequency radiation field (cf. 



|12[ ) . Thus the influence of the LF gravitational radia- 
tion field on the scale factor is small for sufficiently early 
times compared to other forms of radiation. The grav- 
itational wave background considered here is alone not 
able to support a cosmological model: Since the time or 
scale factor dependence of energy density and pressure 
is already fixed, the two equations ([l2]) and ( |l3| ) are not 
compatible. But we may assume that a relativistic fluid 
with a p m — p m /3 equation of state is also present. This 
ensures compatibility at the price of fixing the density 
ratio between fluid and waves and gives 



_ 15a' 2 - 4a 2 fcj 5 07r 

a" + -nkfqa = 0. 
5 

From the second relation the scale factor follows as 
a = aisin(lrj), I = 4irk{q/5, 



(78) 
(79) 

(80) 



thus the energy density of the relativistic fluid changes 
as 

fcfg(3-4sin 2 (/?7)) 



Pm — 



10Ga 2 sin 4 (Zr/) 



(81) 



It appears that the mere existence of a LF gravitational 
wave background could cause a reversal of the expansion, 
however small its energy density may be initially (that is, 
for times n with In << 1) compared to that of the rela- 
tivistic fluid, note that p g /p m = sin 2 (In) / (3 — 4 sin 2 (In)) . 
For In > 0.8861, p g exceeds p m , for nl > it/3 the fluid 
energy density becomes formally negative through the 
interaction with the gravitational wave fields. The aver- 
aged Weyl tensor given by 



(C0101) — (C1212) — 



lQi:a\k\q 
15 



sin 2 (In) 



(82) 



vanishes at the singularity n — and is non-zero for 
n > 0, suggesting that the presence of a LF component 
of gravitational radiation is not a coordinate effect. 

The case of a nonzero limit eo for k — > is more 
complicated and seems not to allow an analytic solution. 
As another example we treat the interaction of the LF 
background with a A— term. It is common practice |25| 1 
to write the A— constant on the right-hand-side of the 
field equations, i.e. to consider a fluid with the com- 
ponents p m = g^g, p m = -gTg. Considering again 
the case lim fe ^ e o = 0, we arrive at ([77|) as before, but 
the field equations ( |l2| ) and ( |l3|) lead now to an expo- 
nential increase or decrease of the scale factor (in the 
77-coordinate) : 



a\ exp (±rjJ Aiik\q/3 



(83) 



As in the former case of a relativistic fluid, the 
"A— matter" and the LF gravitational wave background 
are coupled since 



45GA 



(84) 
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is required for consistency. The Weyl tensor components 

I6na 2 klq 



(C 



0101/ 



(C : 



1212/ 



15 



(85) 



are also finite (and nonzero) at ?/ = 0. 



V. DOMINANT GRAVITATIONAL WAVE 
BACKGROUND 

Let us now assume a pure gravitational radiation uni- 
verse. In the high frequency approximation all is already 
done: We do not have to face an appreciable back re- 
action term in the wave equation, the spectral shape of 
the waves is time-independent, the amplitude decreases 
only due to redshift effects, the equation of state can 
well be approximated by p g — 3p g , and (|12|) , (|13|) (with 
Pm = Pm = 0) give the Tolman radiation cosmos. 

The presence of only LF components in the spectrum 
does not allow a solution based on gravitational waves 
alone, thus the spectrum will not be restricted subse- 
quently. We use the two equations ([30]) and (|5l]) for 
the frequency integrated quantities g^ and f 4 . Adding 
(O),(03) with p m — and p m — 0, one has two further 
equations to determine the four unknown time functions 
/2,/4 7<?4 and a. We solve the last equations for g^ and 
/4 and use the result in (pl|) to find / 2 . This leads to 



g 4 = — (aa" + 3a' 2 ) + — — / 2 - ~-/ 2 , 

87r 6 a z 6 a 



1 / n n m 1° a ' 2 r 
/* - 5-(-<w» + 3a' 2 + ^r— h 
sir 6 a z 
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3 a 
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2- 



(86) 

(87) 
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/^^" + (2(1-2 C )- + -(1-6C)^)/: 
+ 3a /2 ' 



= 2^((3C + 2)^ + (3C-4)^)/ 2 



-,/2 



+ ((C-2)^ + (C + 4)^)/^ 
a a^ 



(89) 



This system may be studied in two cases, (i) assuming 
C = 0, that is, no back reaction in the wave equation and 
(ii) full back reaction, C, = 1. Furthermore one has to 
solve the generalized wave equation ([14]) and check the 
compatibility of its solution with the moments derived 
from (|8q)-(|89|). Finally, if all this succeeds, the energy 
density and pressure of the waves follow from the familiar 
equations 

p g = 3a'7(87rGa 4 ), p g = (-2aa" + a' 2 )/(87rGa 4 ). (90) 



A. C, = 0: Tolman universe 

With C, — 0, ( JS9[ ) can be written as the product of two 
factors: 



(2a'/ 2 - a/a) (aa" - 2a' 2 ) = 0. 



(91) 



Thus two different cases emerge, depending on which fac- 
tor vanishes. If the first factor in (^Tj) is zero, one obtains 
with an integration constant c (which may be gauged to 
zero) 



h = ca 2 



(92) 



This relation is compatible with ( |88| ) if and only if a" = 
0, that is, if the scale factor has the time dependence 
a = brj of the Tolman radiation cosmos J26| , with the 
energy density and pressure given by 



p g = 3p 9 



8irGb 2 r] 4 



(93) 



From (Bq, pTT) it follows that the other moments are time- 
independent: 



.94 = b 2 (c+ —), / 4 



3^ 

8tt' 



(94) 



The average of the Weyl tensor component Cbioi or C1212 
comes out as the constant 2b 2 . The expressions for /a, /j 
and 174 found as solutions of differential equations must 
be compatible with those derived directly from the the 
spectral density ( |70[). The three frequency dependent 



functions l,m,n in ( |70| ) must therefore be chosen so that 
the moments ji,f^-,g\ have the time dependence which 
we have just derived. Self-consistency is only ensured if 
the functions l,m,n exist. Working with real quantities, 
the spectral density / is given by 



/ = 2(n + I cos(2krj) — m sm(2krj)) 



(95) 



in the radiation cosmos. As one verifies from the defi- 
nition of 54 in jigj ) and of f and /1 in (jig), 34 can be 
written 

.94 = 4n 4 - / 4 . (96) 

From the definitions of / 2 and f 4 in (fig ) one obtains 
together with (|f) and (fj) 



k 2 (lcos(2krf) - m sm(2kT]))dk = cb 2 f] 2 /2-n 2 , (97) 
k 4 (lcos(2kr]) - msm(2kr]))dk = -cb 2 /A. (98) 



( pq ) can be obtained from ( |97| ) by differentiation with 
respect to 77, so only ( B7P is needed. All fc-integrations 
considered so far run from fc = to infinity. We formally 
extend l(k) and m(k) to negative values by 



l(k) = l(-k), m{k) = -m(-k) for k < 0. (99) 

This allows us to rewrite ( |97| ) as complex Fourier trans- 
form 



i 



'2nJ- 



(f)(k) exp(2i/c?y) dk 



2^2 



2 c6 ry 

7T 2~ 



"-2J 



(100) 



applied to the complex valued function <fi(k) — k 2 (l(k) + 
im(k)). Since 4>*{k) = <p(—k), the Fourier transform of 
4>(k) is real. To find <f>(k), one has to invoke the Fourier 
inversion theorem. This requires to consider the right- 
hand-side of ( |10C| ) also for negative values of rj. Since 
the right hand side is not absolutely integrable over the 
whole time axis, <p(k) must be understood as generalized 
function (see, e.g., J27), p8| , [ g9j for a confirmation of the 
subsequent calculations). Extending the function space 
in this way, the Fourier inversion theorem remains valid. 
For instance, a polynomial in r\ gives rise to the Dirac 
delta function S(k) and its derivatives in Fourier space. 
One obtains 



k 2 l{k) 



cb 2 d 2 S n 2 
4 dk 2 4-7T 



S(k) 



k m 



(k) = 0. 



(101) 
(102) 



We may also extend n(k) to negative values by n(k) = 
n{—k) for k < 0. The spectral components of energy 
density p(k,rf) and pressure p(k,rf) are then symmetric 
functions of k, and the frequency integrated total density 
can be written p = j/_ dkp(k,rj), with a similar ex- 
tension of the integration interval for the pressure. This 
allows to handle terms involving delta functions apply- 
ing the usual rules for these functions. Replacing second 
derivatives of the delta function using the formulae (cf. 
p7[ , the prime here denotes the derivative with respect 
to k) 

s(k)8"(k) = s"(0)S(k) - 2s'(0)S'(k) + s(0)S"(k), (103) 

one obtains for the spectral decomposition of the energy 
density and pressure from (|52|), (p3l) 



a 4 Gp 9 (£:,?7) = 2n(fc)fc 4 - -^n{k)k 2 + -^-cb 2 5"{k) 



'/ 



(il n2 _ cb 2 )5(k), 

7f 



(104) 



a Gp g — 3a Gp g 



2tlA 



b^c 
~~2 



(106) 



in agreement with (|96|), showing the self-consistency of 
the calculation. 

The singularities of the spectral decomposition p g (k, rf) 
show that the infrared mode k — contributes a finite 
and time-dependent amount Ga A pi r = -^n 2 — ^|- to the 
total energy density. At superhorizon scales, more pre- 
cisely at scales with kr\ < %/3.5, the spectral components 
become negative. Similar conclusions hold for the pres- 
sure, which becomes negative for kr\ < V2.5. The in- 
tegrated values of energy density and pressure however 
stay always positive thanks to the infrared behaviour of 
their spectral components. 



B. ( = 0: de Sitter scale factor 



The vanishing of the second factor in (pll) gives the 
scale factor of the de Sitter universe, a = -4-. The 



density and pressure calculated from (90) are time- 
independent: 



Pg = ~Pg = 



3H 2 

8^G' 



(107) 



The differential equation for f 2l (pq), has the general 
solution 



t r . 7/3 3 ln M 

h = -, + soil 1 - 



r/* 



UttH 2 if 



(108) 



with two constants r$ and sq, where tq may be gauged 
to zero. From (J8qJ87]) one obtains 



, _ 65 1/3 27 1 

•' 4 = ~7T Sor l ~ 7771 tj2 4 - 

9 104 7T_H z! ?7 4 



(109) 
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1° 4. ^W/3 _ ±J^± 4- il^_ (HO) 



'] 



-soV 
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and the Weyl tensor components arc found as 

728 10 1 

(CoiOl) = (C iai2 ) = -^ W /3 ~ ^^4- (HI) 

We have again to check the compatibility of the spectral 
density / (defined by (f70|)) with the time dependence of 
its moments f 2 , f±, g^ as given by the last three equations. 
The integral f 2 corresponding to the definition in (Uq) is 



3a 4 Gp„(fc,77) = 2n(fc)fc 4 - %n{k)k 2 + --\cb 2 6"(k) 

rj z ir] z 

+ {^n 2 -cb 2 )5(k). (105) 

rj z 

The spectral quantities can be integrated immediately 
to give the finite total values 



f 2 = 2n 2 + -^n + 2 k z (l cos(2krj) - m sin(2kr])) dk 



V 2 Jo 



V Jo 



(—1 cos(2krj) + m s'm(2kr))) dk 



k(m cos(2krj) + I sin(2kr))) dk. 



(112) 
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A similar expression holds for fa . Straightforward calcu- 
lation shows that the integrand of 54 = L k 4 g(k, 77) dk 
defined by (f49|), has the form (x = krf) 

g = u cos(2x) + v sin(2.T) + w, (H3) 

2 
u = —r(-l(x 4 - 3.t 2 + 1) + 2mx(x 2 - 1)), 

a; 4 

v = -r(m(x i - 3x 2 + 1) + 2lx(x 2 - 1)), 
x 4 

w = 2n(l - l/x 2 + 1/x 4 ). 

Rewriting the just derived integral expressions for fa and 
#4 as well as for fa as complex Fourier transforms, we 
obtain: 

2 1 2i 

fa = 2n 2 + -^n a + ijj 2 - -ttVo + —ipi, (114) 

77 z ry z 7j 

2 1 2i 

/ 4 = 2n 4 + — n 2 + ip4 5-^2 H "03, (115) 

rj z i~i z 77 

2 2 2i , 

.94 = 2n 4 5-712 + ^-n - ^4 03 

77 z r) 77 

3 2i 1 

H — jV>2 + -o'Vi 700, (116) 

7j z ry 3 77* 

where ipj is a family of time functions defined by 



(117) 



ipj = k 3 (I + im) exp(2ikrj) dk 



(ip2 was already introduced in Section IV) . As in the Tol- 
man case we have extended the functions I, m, n to nega- 
tive values of k by l(—k) = l(k), m(—k) = — m(— k) and 
n(—k) = n(k), thus the functions ipj with even (odd) j 
are real (pure imaginary) . The members of a ^-family are 
connected by differentiation and integration according to 
the rule 



%//j = 2iip 3+1 , 



(118) 



which holds also if the ipj are generalized functions. 
The aim is again to obtain the complex spectral den- 
sity l(k) + im(k) from any of the time functions ipj by 
invertin g the corresponding Fourier integral. We start 
solving (115) for ip 2 , where fa is replaced by the expres- 
sion (109). Using the rule ( 118 ) repeatedly, one obtains 
the differential equation 



r, 2 iP% - 4777^ + 402 



8n 2 + 8n 4 i] 2 + 



27 



260 



267rilV 



-soV 
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(119) 



for ip2- Its solution is given by (note that adding a solu- 
tion of the homogeneous equation would give the wrong 
time dependence) 



1P2 = 2n 2 - 477 4 7j 2 + 13s 7j 7/3 + 



52ttH 2 



V 



(120) 



If one member of a 7/;- family is known, other can be found 
by differentiation and integration. Straightforward cal- 
culation shows, that the T^-functions derived from (12C) 
satisfy also (114) and (116), if fa and 54 on the left-hand- 
sides are substituted from ( |108| ) and (110). We now apply 
the Fourier inversion theorem to (117) with j = 2 and ip 2 
taken from (12C). This requires a continuation of the 
(real) function ip 2 into the region 77 < 0. Here only the 
term proportional to sq requires more consideration, but 
this term contributes nothing to p g and p g , we can there- 
fore put the integration constant sq equal to zero. One 
then obtains 



k l(k) = 77,4 
k 2 m{k) = 0. 



d 2 5(k) 
dk 2 



2n 2 S(k) - 



m 

26ttH 2 ' 



(121) 
(122) 



Again infrared modes enter the spectral density l(k), but 
contrary to the Tolman case one is free to specify the 
spectral function n(k). If n(k) is chosen as zero, the delta 
function terms in / and in the spectral decomposition of 
p and p vanish, but / has still singular terms seen in an 
expansion around k = 0: 



/ 



1 



27/ 4 

13nH 2y k ' rfk*' 39nH 2 



r<7 + 



k 3 + o(k 5 ) (123) 



(we missed in section IVC the de Sitter case, since we 
had excluded infrared singularities). In spite of this sin- 
gularity, density and pressure of the gravitational waves 
integrate to the finite constant values ( |107| ). Additional 
terms from a spectral function n(k) 7^ add nothing 
to the total values p g and p gi thus the singular infrared 
(k = 0) component in p g is cancelled by the integrated 
contribution of 77-modes with k ^ 0, the same holds for 
the pressure. 

It is easy to extend the calculation to account for 
a cosmological constant by adding matter terms p m — 
A/(87rG) and p m = —p m to the equations ([l^Jl3] ). Re- 
peating the calculation at the beginning of this section, 
one obtains (for Q — 0) the same product (|9l]) as in the 
absence of a A constant. The vanishing of the first fac- 
tor gives a Tolman-de Sitter model, where gravitational 
radiation has an equation of state (EOS) of the form 
Pg = Pg/3 an d an energy density decaying as a -4 . The 
scale factor 0,(77) can be expressed in terms of an elliptical 
integral. 

More interesting is the model corresponding to the van- 
ishing second factor, since here the de Sitter scale factor 
follows. No new calculation is needed: In equation (107) 
we have to substitute for p g and p g the total values p g +p\ 
and p g + p\. In the calculations following this equation 
we must only replace the de Sitter Hubble constant H 
in all equations by H = H/{\ - 64ir 2 G 2 A/(3H 2 )) 1 / 2 . 
The de Sitter expansion is generated by two independent 
"matter" sources, by a genuine A constant as well as by 
a suitable spectrum of gravitational waves. We may have 
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an arbitrary (but time-independent) mixture of both in- 
gredients. For sufficiently small A gravitons dominate. 
If A reaches the threshold A* = H 2 /(64tt 2 G 2 ), p g be- 
comes zero and turns to negative values for still larger 
A, to ensure the same total energy density for a different 
composition. 

We have not discussed in this article the origin of the 
primordial wave spectrum, but it should be noted that 
the expressions for energy density and pressure intro- 
duced here may be of interest for concrete models. In 
the case of a quantum origin due to vacuum fluctuations 
in a de Sitter cosmos Q, the produced gravitons can 
be described by a two-point correlation function, which 
for Bunch-Davies vacuum |3l]] corresponds to a classical 
spectral density / given by 



/, 



BD 
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kir 2 rfk 2 
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(124) 



Ibd is obtained by comparing the quantum expectation 
values of bilinear terms in the metric (as given, e.g., in 
|p2fl ) with the stochastic averages discussed here. The 
spectral densit y Jbd has the same infrared singularity as 
the expression (|l23| ) derived previously. Since the coinci- 
dence of both spectra holds only approximately for small 
k, the back reaction of the Bunch-Davies gravitons on the 
scale factor will change the scale factor^]. Again, in spite 
of the singularity of the spectral density /bd, the inte- 
grated values of energy density and pressure require no 
infrared cut-off. If one introduces an ultraviolet cut-off 
at the frequency k\ with x± = kirj, one obtains for energy 
density and pressure of the Bunch-Davies gravitons 



Pa 



HH 4 



x\(x\ - 7), p g 



hH 4 



xfixf + l). (125) 



Considering limiting cases, for high frequencies ii > 1 
follows the expected EOS p = p/3, for low frequencies 
iiCl one has p = —p/i (together with p < 0). The 
expressions (125) for p g and p g apply only for a de Sitter 
scale factor. As a result of to back reaction, the local 
values of p g and p g change with the background geome- 
try (thus they are, in a sense, not local): The lower the 
frequency, the stronger is the dependence of the effective 
EOS on the background gravitational field. 



C. f = 1: Full back reaction models 

Unfortunately, the approach described in the previous 
subsections does not work in the real case £ = 1. Even 



^Similar conclusions follow from the important work by 
Tsamis and Woodard, who have discussed quantum gravity 
back reaction on an inflationary expansion rate in numerous 
papers 0. 



the first step, finding analytic solutions of (|88|) , (|8£|) for 
a, /2, has been so far unsuccessful. The back reaction 
function can be written b = C(ir + ^~) ^ or a domi- 
nant gravitational wave background, thus the differential 
equation for the wave amplitudes h takes for C = 1 the 
unusual form 



n" n' 2 

h" + h(k 2 + — + 2%) = Q. 



(126) 



This equation suggests that graviton creation persists in 
a wave dominated universe, provided the scale factor is 
different from a ~ 7/ 1 / 3 . Note that the latter scale factor 
is inconsistent with the system (|88|) , (p9|) . For a treat- 
ment of the general case one has to resort to numerical 
calculations, which will be discussed elsewhere. 



VI. FINAL REMARKS 

We have seen that the stochastic back reaction equa- 
tions treated here form, on the one hand, an apparently 
self-consistent system of equations with interesting solu- 
tions. On the other hand, it is not clear, how far the 
solutions deviate from true solutions of Einstein's field 
equations, either for some range of parameters or in some 
regions of space-time. We shortly discuss what could be 
done to clarify and to improve the situation. 

One has to realize that ensemble averages are consid- 
ered, which always differ from actual realizations of a ran- 
dom process. This is inherent to the method and cannot 
be changed, but one is able to say something more about 
statistical deviations from true solutions, if a Gaussian 
or some other process is assumed. 

The main shortcoming of the approach is the use of 
a second-order approximation to general relativity, but 
improvements are possible. In principle, Monin and Ya- 
glom's statistical treatment of nonlinear field theories 
works for arbitrary nonlinearities, if they are present in 
polynomial form. Only the technical complexity grows 
in a full treatment, since many higher-order correlation 
functions must be taken into account. Writing the Ein- 
stein field equations as 



1 



(R^ ~ ^uR pa 9 P <r)g a = KT^g* 



(127) 



where g is the determinant of the metric tensor, we see 
that the left-hand sides consist of huge polynomials of 
maximal degree 12 in the (exclusively) covariant compo- 
nents of the metric tensor and its first and second deriva- 
tives. Taking the ensemble average of these expressions 
leads to correlation functions up to sixth order. If the 
random process g^ v is Gaussian, the higher order corre- 
lation functions can be reduced to the second-order func- 
tions studied in this article. This would allow us to turn 
the back reaction equations into -in some sense - exact 



12 



relations. Present discussions on a primordial stochastic 
gravitational wave background usually assume a quan- 
tum origin of gravitons, which are born out of zero-point 
vacuum fluctuations. The classical correlations discussed 
here are expected to be related to quantum mechanical 
expectation values, thus it seems natural to assume Gaus- 
sianity. 
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